Twisted bilayer WSe\(_2\) DFT

Lattice structure for \(\sim 22\) degrees

We perform a DFT calculation of twisted bilayer WSe\(_2\) at a twist angle of approximately 22 degrees. For a twist angle there are two possible ways of stacking the layers due to lack of \(C_6\) symmetry of the monolayer:

  • Parallel stacking (P-WSe\(_2\)): here the layers can be translated in such a way that two metal atoms and two chalcogen atoms are vertically aligned. Inversion symmetry is broken for this stacking.

../../_images/logo.png

Fig. 1 P-WSe\(_2\) where the cell border vertices correspond with two aligned metal atoms. Two chalcogen atoms align within the circle of metal atoms.

  • Anti-parallel stacking (AP-WSe\(_2\)): here no two chalcogen atoms align and inversion symmetry is obtained, a symmetry the normal monolayer does not have.

../../_images/logo.png

Fig. 2 AP-WSe\(_2\) where the cell border vertices correspond with two aligned metal atoms. No two chalcogen atoms align perfectly but each circle of metal atoms contains one chalcogen atom, allowing for inversion symmetry.

DFT

For the dft calculation we use the P-WSe\(_2\) configuration of the bilayer. Before performing the calculation of the bands we perform a relax calculation. For both we use the software package Quantum Espresso (QE). QE also offers the possibility of projecting the standard atomic orbitals onto the wavefunctions. This we do to gain insight of the nature of the highest valence bands, which form the flat band at the magic twist angle.

../../_images/logo.png

Fig. 3 The band structure of P-WSe_\(2\) for a twist angle of \(\sim\)22 degrees. Like for the monolayer spin-orbit coupling (SOC) splits the valence band at the \(K\) point. We show the contribution of the orbitals where we discern for the layer and spin nature.

In figure Fig. 3 we plot the bands with the percentile contribution of the atomic orbitals. The most striking feature of these bands is the degenerate nature of the top valence band. Degeneracy indicates that we can recombine the wavefunctions of these two bands such that each lives on its own layer and does not couple with the other. If it is this band that constitutes the flat band for smaller twist angle we can conclude that interlayer hopping does not contribute to the forming of the flat band.

The valence bands at the \(K\) point with lower energy however do show clear signs of interlayer coupling. These bands are non-degenerate and as such must live on both layers. It is however unlikely that backfolding or interlayer coupling can compensate for the 0.5 eV deficit required for these bands to end up on top.